Inhomogeneity— induced enhancement of the pairing interaction in cuprates 
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Scanning tunneling spectroscopy has recently discovered a positive correlation between the mag- 
nitude of the superconducting gap and positions of dopant oxygen atoms in Bi-based cuprates. 
We propose a microscopic mechanism that could be responsible for this effect. In particular, we 
demonstrate that the dopant-induced spatial variation of the atomic levels always enhances the 
superexchange interaction. 
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Scanning tunneling microscopy (STM) and scanning 
tunneling spectroscopy (STS) have recently confirmed 
that nanoscale electronic inhomogeneity is an inherent 
feature of many groups of high-temperature supercon- 
ductors (HTSC). Being capable of direct probing the 
local density of states (LDOS), these methods revealed 
strong spatial modulation of the energy gap in Bi-based 
compounds This modulation occurs on a very short 
length scale, of the order of the coherence length. It 
has been shown that the system consists of regions of 
relatively small gap (A w 25 ~ 35meV) with high and 
sharp coherence peaks and regions of larger gap (A w 
50 -r 75meV) with small and broad peaks [2|. Since its 
discovery, the inhomogeneity was commonly attributed 
to a disorder introduced by poorly screened electrostatic 
potential of the out-of-plane oxygen dopant atoms. A lot 
of theoretical works exploited this idea, predicting modi- 
fication of the energy gap in the vicinity of a dopant atom 
Q ■ Very recently, STM experiments have shown a strong 
correlation between position of the dopant atoms and all 
manifestations of the nanoscale electronic disorder 0, [B| . 
Thus, these experiments proved the impurities to be the 
source of the inhomogeneity. On the other hand, they 
revealed a very important feature: there is a positive 
correlation between the position of a dopant atom and 
the magnitude of a gap [4| . The sign of this correlation 
function contradicts the previous theoretical predictions, 
based on a direct reduction of the gap by the dopant's 
electrostatic potential Q. Therefore, a new mechanism 
that is capable of gap enhancing close to impurities, is 
needed to explain the correlations. It was shown by Nun- 
ner et al. |6| , that an assumption of an enhancement of 
the pairing potential by the presence of a nearby dopant 
atom leads to the correct correlation between the height 
of the coherence peaks and the magnitude of the gap as 
well as between the position of the dopant atom and the 
magnitude of the gap. The spectral properties could also 
be explained in a different way, assuming that the ob- 
served peaks at the edges of the gap arise from resonant 
bound states rather than they are the coherence peaks 
0- 

From the theoretical point of view, various parame- 
ters of the microscopic models of HTSC can be modified 
by the presence of the dopant atoms. These atoms are 



charged impurities and since cuprates are close to the 
insulating state, the number of carriers is too low to ef- 
fectively screen their electrostatic potential. As a result, 
atomic levels in the Cu02 plane are shifted in the vicin- 
ity of the dopants. The presence of the dopants also 
induces a local distortion of the lattice 0. This, in turn, 
may modify the hopping matrix elements as well as the 
electron-phonon coupling. It has recently been shown 
that inhomogeneity itself can increase the superconduct- 
ing transition temperature [^, [l(| • 

Here, we show that in strongly correlated systems, the 
position-dependent shift of the atomic levels alone is suf- 
ficient to enhance the pairing interaction, thereby leading 
to the correct sign of the dopant-gap correlation func- 
tion. We assume that HTSC can be described by the 
t-J model, with the exchange interaction as the main 
pairing mechanism . The physical picture behind our 
analysis is based on the same arguments, which allow 
one to derive the t-J model from the Hubbard one 
Since strong Coulomb repulsion excludes double occu- 
pancy, eigenstates of the half-filled Hubbard model are 
built out of localized electron states, with one electron 
per lattice site. A virtual hopping of an electron to a 
neighboring site (occupied by an electron with the oppo- 
site spin) gives rise to an effective spin exchange inter- 
action with potential J = 4t 2 /U, where t is the hopping 
integral and U is the Coulomb potential. Obviously, local 
modification of the hop ping matrix elements leads either 
to an enhancement [1 31 ] or to a reduction of the exchange 
interaction. Recent density functional investigations in- 
dicate that the dopant oxygen atoms are responsible for 
displacements of atoms in the Cu02 plane [|| . Since these 
displacements are perpendicular to the plane, such a dis- 
tortion increases the interatomic distances. Therefore, 
one might expect a reduction of the hopping integrals 
and a simultaneous decrease of J. However, the effective 
exchange interaction should depend also on the atomic 
levels at sites, which are involved in the virtual hopping. 
It comes from the fact, that the energy of the virtual 
doubly occupied state depends on these levels. The basic 
aim of the present manuscript is to show that the diag- 
onal disorder always leads to an enhancement of the an- 
tiferromagnetic exchange interaction. As a result the su- 
perconducting gap increases in the vicinity of the dopant 
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atoms. 

We start with the two-dimensional one band Hubbard 
model on a square lattice 

H = -tJ2 ( a L a 3° + h - c ) + - m)4t°^ 

(i.j)a- icr 

a i1; a iT a il a ilj (1) 

i 

where a\ a creates an electron with spin a at site i, /i 
is the chemical potential and Vi is the atomic level at 
site i. The latter quantity accounts for the electrostatic 
potential of the dopant atoms. 

In the following we argue that the enhancement of the 
pairing interaction is driven by a random distribution 
of Vt. However, there is a severe drawback of such as- 
sumption: strong enough dopant's potential should sig- 
nificantly affects the electron concentration in its neigh- 
borhood, resulting in a highly inhomogeneous charge dis- 
tribution. On the other hand, the STS experiments show 
that the root mean square of the electron density in 
these systems are below 10% 14]. One can expect how- 
ever, that strong Coulomb repulsion along with vicinity 
to half-filling, meaningfully reduce the dopant-induced 
inhomogeneity. Namely, the decrease of the local den- 
sity of electrons close to the dopant atoms must result in 
its increase away from them. Such an increase is ener- 
getically highly unfavorable. In particular, at half-filling 
and for a strong Coulomb repulsion impurities cannot 
produce any charge inhomogeneities at all. In order to 
verify the effectiveness of this mechanism for a physi- 
cally relevant doping, we have performed exact diagonal- 
ization study on a finite-size cluster. We have consid- 
ered two impurities located z = 1.5a above the cluster 
with a being the lattice constant. It is the distance be- 
tween the dopant oxygen atom and the Cu02 plane in 
Bi-based HTSC. These impurities act as a source of a 
screened electrostatic potential, that shifts the atomic 
levels Vi = VbX] m ex P(^^mi/A)/fl mi , where the sum- 
mation is carried out over all dopant atoms. Here, R m i 
is a distance between site i and site m (above which the 
dopant atom is located) and R mi = s/R^i + z 2 . Such 
a form of Vi accounts for the fact that the potential is 
screened only by electrons moving in the Cu02 plane. 
However, our final conclusions are independent of the 
form of V, . 

Fig. Q] presents root mean square of the local electron 
concentration 5n as a function of the Coulomb interac- 
tion U for various values of Vq. One can note a signifi- 
cant reduction of the charge inhomogeneity in a large-t/ 
regime in spite of the strong variation of the electrostatic 
potential V,. In the uncorrelated case, inhomogeneity of 
the charge distribution depends on the ratio Vo/t. Con- 
trary to this, in the strongly correlated nearly half-filled 
system inhomogeneity is determined predominantly by 
Vq/U . Therefore, small charge density variation observed 
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FIG. 1: Root mean square of the local electron density as a 
function of U for various values of Vo, as explicitly indicated 
in the figure. Results have been obtained for 10 electrons on 
a 12-site cluster with periodic boundary conditions. A equals 
two lattice constants. 



in STS experiments does not require weak electrostatic 
potential of the dopant atoms. 

We now turn to the main problem, i.e., the enhance- 
ment of the pairing interaction by the dopant atoms. 
For this sake, we make use of the Hubbard operators 
X^' 13 = \a)i(j3\i {a, (3 = 0,T,j,U) and carry out the 
canonical transformation H — > exp(— S)H exp(S'). The 
generator S = —S^ is chosen in such a way that the 
transformation eliminates the transfer of electrons be- 
tween lower and upper Hubbard subbands. Straightfor- 
ward calculations lead to the generator of the form 



S= E s ^ 



X l ' a X 



U + (Vi-Vj) U-(Vi-Vj) 



, (2) 



where s(f) = 1 and s(l) = — 1. 

Next, one can project out states with doubly occupied 
sites and obtain the effective Hamiltonian acting in the 
space spanned by states |0)i, | |)< and | J,)*. This is a 
generalized t—J Hamiltonian with a site dependent spin 
exchange interaction: 

Ht-J = -t (oLojv + h.c.) + ^2(Vi - li)a\ a a i<T 



^ Jij [Si ' Sj jn-iUj 



(3) 



Here, a ia creates an electron at site i if this site previously 
had no electron. The interaction Jij is given by 
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(4) 
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where 



Vij 



{Vj-Vjf 



u 2 - (yi - v,) 2 



> 0. 



(5) 



In the homogeneous case (Vi = const), J t j reduces it- 
self to the standard form J = At 2 /U . Otherwise, it 
is always larger than Jo, independently of the distribu- 
tion of the atomic levels and any particular form of Vi. 
Moreover, the strength of the coupling increases with the 
difference between the atomic levels at neighboring sites. 
Assuming Jij to be the effective pairing interaction one 
comes to the conclusion that superconductivity can be 
enhanced by inhomogeneous distribution of the atomic 
levels. Since these levels directly affect the pairing in- 
teraction, spatial variation of the superconducting gap 
does not require strongly inhomogeneous charge distri- 
bution. The surprising enhancement can be easily un- 
derstood from the analysis of virtual processes presented 
in Fig. [2l The resulting exchange interaction is propor- 
tional to the squared hopping matrix elements (t 2 ) and 
inversely proportional to the energy of the intermediate 
state with doubly occupied site. There are two second or- 
der processes that lead to the effective antiferromagnetic 
coupling of spins at sites i and j. The first of them (Fig. 
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FIG. 2: (Color online) Second order two-site processes con- 
tributing the effective coupling between spins at sites i and 
j- 

Hi) gives J° = 2t 2 /[U + (Vi - Vj)], whereas the second 
one (Fig. [2p) gives = 2t 2 /[U-{Vi - Vj)]. Taking into 
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account both these processes one obtains 
what is the value given by Eq. [3J 

Since the Coulomb repulsion is the dominating energy 
scale in cuprates, one can assume that \Vi — Vj\ <C U. 
Then, an approximate formula rjij ~ (Vi — Vj) 2 /U 2 holds 
and one can see that r\ij decreases much faster than Vi 
when the distance from the dopant atom increases. As 
a result, a significant enhancement of the exchange in- 
teraction is expected only in the closest vicinity of the 
dopant. 

It is known that an enhancement of the pairing poten- 
tial close to impurities leads to the correct sign of the 
impurity-gap correlation function as it has been shown 
in Refs. [6|, [l4|. However, in these papers the spatial 



origin. Contrary to this, we have proposed a microscopic 
mechanism that gives an explicit form of this dependence. 
In the following, we investigate whether the proposed sce- 
nario quantitatively reproduces the STS data. In partic- 
ular, it is necessary to check whether a significant en- 
hancement of the superconducting gap can be obtained 
for the model parameters, which do not lead to strong 
fluctuations of the electron concentration. 

In order to estimate the spatial variation of the super- 
conducting order parameter, we have assumed the res- 
onating valence bond (RVB) scenario ll| and investi- 
gated the Hamiltonian ((3|) in the mean-field approxima- 
tion with renormalized hopping integral £—►£=(! — n)t 
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15j . Within such an approach the constraint of no 
double occupancy is not fulfilled exactly. On the other 
hand, it is a method that allows one to investigate super- 
conductivity in inhomogeneous systems by means of the 
Bogoliubov-de Gennes (BdG) equations. Following the 
experimental results we restrict further study to such val- 
ues of the model parameters U and Vb, which give small 
electron density fluctuations Sn < 0.1 (see Fig. [TJ. There- 
fore, we neglect the term containing Vi in the Hamilto- 
nian ([3|) , that is responsible for the charge inhomogeneity. 
Since the mean-field RVB approach leads to an incorrect 
doping dependence of the critical temperature, this term 
could be responsible for an unphysical contribution to 
the spatial dependence of the order parameter. 
The mean-field Hamiltonian takes on the form: 



H 



t-j 



■ E ( a L a 3° + h ' C -) _ ^E a ta a ™ 
(i,j)a io 

E [ A ^- (°Mi + a )Ai) + h - c ] . ( 6 ) 



where A,- 



-Jij(a.iiaj^ + aj^ajf)/2. It can be diago- 



nalized with the help of the transformation: 

Cia — ^ ^ ( M m7rt(T — S ( IJ ) V initio ) ■> 
n 

where Ui n and Vi n fulfill the BdG equations: 

n 



E 



(7) 



(8) 



with TLi 



-ts 



/iSij. Here, S 



1 for the neigh- 



boring sites i,j and otherwise. The superconducting 
order parameter is determined self-consistently by: 



A« = ^E 



3 " 



+ Uj n v?„) tanh 



2k B T 



(9) 



variation of the pairing potential is of phenomenological 



We have solved the BdG equations for the d-wave sym- 
metry. The calculations have been carried out on a 32 x 32 
cluster with 170 randomly distributed impurities, what 
gives concentration 16%. Assuming that each dopant 
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FIG. 3: (Color online) A t /t calculated for U = 8t, V = 5t, 
X — 2a and n ~ 0.83 The positions of impurities are marked 
by dots. The dashed line shows the cross-section, along which 
the LDOS presented in Fig. [4] was calculated. 




FIG. 4: (Color online) LDOS calculated along the line marked 
in Fig. 0for the same model parameters. 



neighborhood of an impurity. Consequently, the pair- 
ing interaction is enhanced in a region larger than the 
coherence length only provided this region contains sev- 
eral impurities. Figs. [3] and |4] show that the mean-field 
approach correctly reproduces also the main qualitative 
features of the local density of states (LDOS). In partic- 
ular, increase of the superconducting gap in the vicinity 
of impurities is accompanied by a reduction of the height 
of the coherence peaks. The reduction is however much 
weaker than observed experimentally. This discrepancy 
can be attributed either to the limitations of the mean- 
field analysis or/and to fact, that we have considered only 
the increase of the pairing interaction and neglected other 
possible dopant-induced effects, e.g., modulation of the 
hopping integral. 

To conclude, we have derived a purely microscopic 
mechanism that can be responsible for the observed en- 
hancement of the superconducting gap in the vicinity of 
impurities. Our approach follows the generally accepted 
view that high temperature superconductivity originates 
from the exchange interaction between electrons in the 
lower Hubbard subband. We have considered a simple 
scenario of the influence of dopant atoms on electrons 
moving in the Cu02 planes. It is based on the shift of 
the atomic levels by the dopant's electrostatic potential. 
One of the most attractive features of this approach con- 
cerns its simplicity. On the other, hand such a simple 
approach neglects other possible effects originating from 
the presence of the dopant oxygen atoms. In particular, 
we expect that vibrations of these atoms can further re- 
duce the height of the coherence peaks in the large gap 
areas, as observed experimentally. 

This work has been supported by the Polish Ministry of 
Education and Science under Grants No. 1 P03B 071 30. 



introduces one hole into the Cu02 plane, we have per- 
formed calculations for the same hole concentration. Fig. 
[3] shows the spatial variation of the superconducting or- 
der parameter obtained for U = 8t, Vo = 5t and A = 2a. 
For such model parameters exact diagonalization study 
gives the root mean square of the electron density below 
0.1. (see Fig. [IJ, what is comparable to the experimen- 
tally observed fluctuations One can note a signifi- 
cant spatial variation of the superconducting gap. We 
have found that the root mean square of the order pa- 
rameter equals 25% of its average value. As expected Aj 
increases in the vicinity of each impurity. However, the 
degree of this enhancement strongly depends on a con- 
figuration of impurities. It is significant in the regions, 
where several impurities are located close to each other. 
Otherwise it is much weaker. This result follows from the 
specific form of exchange interaction as given by Eq. f4|. 
Due to the screening of the electrostatic potential, rjij 
gives nonnegligible contribution to Jy only in the closest 
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